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Introduction: Calculation of oxygen 

fugacity in high pressure and temperature 
experiments in metal-silicate systems is 
usually approximated by the ratio of Fe in the 
metal and FeO in the silicate melt: 
AIW=2*log(X Fe /X Fe0 ), where IW is the iron- 
wustite reference oxygen buffer. Although this 
is a quick and easy calculation to make, it has 
been applied to a huge variety of metallic (Fe- 
Ni-S-C-O-Si systems) and silicate liquids 
(Si0 2 , AI2O3, Ti0 2 , FeO, MgO, CaO, Na 2 0, 
K 2 0 systems). This approach has surely led to 
values that have little meaning, yet are applied 
with great confidence, for example, to a 
terrestrial mantle at "IW-2". Although f0 2 can 
be circumvented in some cases by 
consideration of Fe-M distribution coefficient, 
these do not eliminate the effects of alloy or 
silicate liquid compositional variation, or the 
specific chemical effects of S in the silicate 
liquid, for example. In order to address the 
issue of what the actual value of f0 2 is in any 
given experiment, we have calculated f0 2 from 
the equilibria 2Fe (metal) + Si0 2 (liq) + 0 2 = 
Fe 2 Si0 4 (liq). 

Calculations: Modeling this equilibria at 

high PT conditions requires satisfactory a-x 
relations and thermodynamic or equation of 
state (EOS) data. We utilize activities of Fe in 
the FeNi alloy [1,2] and in the FeNiS system 
[3], and activities of Si0 2 and Fe 2 Si0 4 in the 
silicate melt [4]. EOS data for FeNi alloy [5] up 
to 30 GPa are used, as is the EOS for silicate 
melt [4]. The position of the IW buffer also 
changes with P and T, and is calculated from 
recent work [6]. 

Calculations along hypothetical adiabats 
and for specific experiments: To illustrate 
the potential differences between a basaltic 
melt and a peridotite melt, we can first 
calculate f0 2 for two hypothetical melt 
compositions along an adiabat (0.3 K/km 
gradient) between 1 bar and 50 GPa. 
Calculations have been carried out for a few 
specific experiments from the literature chosen 
to represent certain extremes of composition or 
PT conditions. For example, high T runs from 
[7] and [8] are chosen to investigate the effect 
of very high T (2873-3000 K at 25-26 GPa). 
Sulfur-rich melts of [9] and [10] were chosen at 
very low and very high pressures (0.8 vs. 25 
GPa). And several high PT experiments from 


[11,12] were chosen for additional 
comparisons. 

Results: We can compare the actual 
calculated values to those calculated using the 
ratio approximation approach. In general, the 
calculated values are 1 to 2 log f0 2 units 
higher than those using the ratio 
approximation. The difference becomes higher 
at high pressures (>25 GPa) for the peridotite 
(Fig. 1). The difference stays relatively low for 
the basalt, but the EOS for basalt is more 
uncertain at higher pressures due to the role of 
alkalis. These results suggest that calculated 
f0 2 from the ratio approximation could be too 
low, and the offset becomes even higher at the 
high PT conditions estimated for core 
formation. As a result this should be 
considered in detailed modeling efforts for a 
complete understanding of metal-silicate 
equilibria. 



Fig.l: Calculated AlW (at P and T) vs. pressure for 
peridotite and basalt (blue curves), and individual 
experiments from the references shown in bracketed 
numbers. Numbers in parentheses are AlW values 
calculated using the ratio approach. Experiments of [9] 
and [10] involve FeNiS melts and the difference between 
the ratio approach and new calculations can be close to 5 
log f0 2 units. Similarly, the high PT experiments of [11] 
utilizing MgO-rich melts and FeNi metal are 1-2 logf0 2 
units different. Values for other experiments at higher 
temperatures [7,8], or low pressures [12], are closer, but 
still - 1 logf0 2 unit different. 
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